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Abstract 

It is known that the diffuse Ha emission outside of bright H II regions not only are very extended, but also 
can occur in distinct patches or filaments far from H II regions, and the line ratios of [S II] A6716/Ha and [N II] 
A6583/Ha observed far from bright H II regions are generally higher than those in the H II regions. These 
observations have been regarded as evidence against the dust-scattering origin of the diffuse Ha emission 
(including other optical lines), and the effect of dust scattering has been neglected in studies on the diffuse Ha 
emission. In this paper, we reexamine the arguments against dust scattering and find that the dust-scattering 
origin of the diffuse Ha emission cannot be ruled out. As opposed to the previous contention, the expected dust- 
scattered Ha halos surrounding H II regions are, in fact, in good agreement with the observed Ha morphology. 
We calculate an extensive set of photoionization models by varying elemental abundances, ionizing stellar 
types, and dumpiness of the interstellar medium (ISM) and find that the observed line ratios of [S II]/Ha, 
[N II]/Ha, and He I A5876/Ha in the diffuse ISM accord well with the dust-scattered halos around H II regions, 
which are photoionized by late O- and/or early B-type stars. We also demonstrate that the Ha absorption 
feature in the underlying continuum from the dust-scattered starlight ("diffuse galactic light") and unresolved 
stars is able to substantially increase the [S II]/Ha and [N II]/Ha line ratios in the diffuse ISM. 

Subject headings: Galaxy: halo — H II regions — ISM: structure — radiative transfer — scattering 



1. INTRODUCTION 

The diffuse Ha emission outside of bright H II regions is 
ubiquitous in late-type galaxies and is generally believed to 
probe the warm ionized medium (WIM; also called diffuse 
ionized gas, DIG), which may be a major component of the 
interstellar medium (ISM) of our Galaxy and other late-type 
galaxies (Reynolds 1991; Walterbos & Braun 1996; Dettmar 
2000; Rand 2000; Hidalgo-Gamez 2005; Haffner et al. 2009). 
It has been argued that the WIM is mainly photoionized 
by ionizing radiation (Lyman continuum; Lyc) that leaks 
out of bright H II regions associated with O stars (Mathis 
1986a, 2000; Domgorgen & Mathis 1994; Sembachetal. 
2000; Wood & Mathis 2004; Woodetal. 2010), although it 
is not clear how the Lyc can penetrate the diffuse H I that is 
observed everywhere in the galaxies (e.g., Seon 2009). 

Previous studies of the diffuse Ha and other optical 
emission lines concluded that the effect of dust scatter- 
ing is negligible (Reynolds 1985; Walterbos & Braun 1994; 
Ferguson et al. 1996a,b; Wood & Reynold 1999). However, it 
was recently revealed that the Ha excess intensity in a num- 
ber of high-latitude clouds is due to scattering by interstel- 
lar dust of Ha photons originating elsewhere in the Galaxy 
(Mattila et al. 2007; Lehtinen et al. 2010; Witt et al. 2010). 
More recently, Seon et al. (2011a,b) presented a strong con- 
nection between the diffuse Ha and far-ultraviolet (FUV) con- 
tinuum backgrounds, which likely suggests a common radia- 
tive transfer mechanism responsible for a substantial fraction 
of both Ha line and FUV continuum photons. Through a 
comparison of the Ha and FUV continuum backgrounds, we 
proposed that the diffuse Ha emission at high latitudes may 
be dominated by late O- and/or early B-type stars (hereafter, 
late OB stars) and dust scattering (Seon et al. 201 lb). Further- 
more, Brandt & Draine (2012) detected emission lines Ha, 
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H/3, [N II] A6583, and [S II] A6716 at high Galactic latitudes, 
whose intensities are correlated with the 100 /im intensity, 
providing an independent confirmation of the importance of 
dust-scattered components in the optical emission lines. If the 
dust scattering indeed significantly contributes to the diffuse 
Ha emission and the effect of dust scattering is not properly 
taken into account, the energy requirement for the ionization 
and volume filling fraction of the WIM will be overestimated. 
Also, the widely-applied practice of using the all-sky map of 
diffuse Ha intensity (Finkbeiner 2003) as a template for both 
the intensity and small-scale structure of the Galactic fore- 
ground free-free emission in cosmic microwave studies (e.g., 
Dickinson et al. 2009; Gold et al. 201 1) needs to be reevalu- 
ated, if dust scattering contributes significantly to the intensity 
and spatial structure of diffuse Ha light outside the actual H II 
regions. 

The recent findings about the importance of dust scatter- 
ing motivated us to reexamine the rationales behind previous 
conclusions against the dust-scattering origin of the diffuse 
Ha emission. It was also necessary to confirm our previous 
results that nebulae photoionized by late OB stars are capable 
of reproducing the typical values of the [N II] A6583/Ha and 
[S II] A6716/Ha line ratios. 

The rejection of dust scattering as a significant contrib- 
utor to the diffuse Ha and other optical line emissions is 
mostly based on two arguments. The first argument is based 
on the Ha morphology in the surroundings of H II regions. 
Ferguson et al. (1996a,b) assumed that in R-bmd continuum 
images, bright OB associations would show /?-band contin- 
uum halos as prominent as the Ha halos around H II regions, 
if the Ha halos around H II regions are caused by dust scat- 
tering. However, they found no corresponding halos around 
bright OB associations in the /?-band images. Additionally, 
light scattered by dust was assumed to be usually concen- 
trated around H II regions, but the diffuse Ha emission is 
often observed in distinct patches or filaments far from H II 
regions (e.g., Hoopes et al. 1996). Consequently, the ex- 
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tended Ha emission outside of bright H II regions was re- 
garded as evidence of leakage of Lyc photons rather than as 
the dust scattering of Ha photons. Second, the optical line 
ratios of [N II]/Ha and [S II]/Ha observed outside of H II re- 
gions have been found to be generally higher than the corre- 
sponding ratios in bright H II regions (Reynolds 1985, 1987; 
Walterbos & Braun 1994). It was thereby concluded that the 
dust-scattering effect could be ignored in spite of the fact that 
bright H II regions are usually surrounded by a large amount 
of molecular gas and dust. 

In the present study, by using radiative transfer models, we 
reexamine whether these arguments are indeed valid in re- 
jecting the significance of dust scattering in the diffuse op- 
tical emission. First, we argue that the assertion based on 
the halo extent or morphology is clearly wrong. Instead, we 
show that the global Ha morphology in face-on galaxies can 
be readily explained by dust scattering. Second, we find that 
the observed line ratios of [N II]/Ha and [S II]/Ha outside of 
H II regions can be explained reasonably well by the dust- 
scattered halos around the H II regions photoionized by late 
OB stars, which is consistent with the main conclusion of 
Seon et al. (2011b). We demonstrate these results by calcu- 
lating photoionization models of the radiation-bounded H II 
regions and putting the resulting emissivity distributions of 
the optical lines into the radiative transfer models as input 
sources for a dust scattering simulation, thereby showing that 
the dust-scattering origin of the diffuse Ha emission cannot 
be ruled out. Lastly, we also point out that the Ha absorption 
line in the dust-scattered starlight (so called "diffuse galactic 
light," DGL) and the stellar continuum of unresolved stars can 
significantly increase the observed values of the line ratios. 

This paper is organized as follows. In Section 2, we dis- 
cuss the morphology and the surface brightness distribution 
of Ha halos around H II regions caused by dust scattering. 
Section 3 investigates one-dimensional and three-dimensional 
photoionization models, and the line ratios resulting from the 
models. In Section 4, we discuss the effect of the underly- 
ing Balmer absorption line on the optical line ratios. Addi- 
tional evidence supporting the results are described in Section 
5. Section 6 summarizes the results. 

2. DUST-SCATTERED Ho HALOS 
2.1. Model 

We performed Monte Carlo radiative transfer simulations 
of the dust-scattered Ha halos around extended sources. For 
the dust-scattering simulation, we used a three-dimensional 
Monte Carlo radiative transfer code (Lee et al. 2008; Seon 
2009), which is similar to those described in Wood & Reynold 
(1999) and Gordon et al. (2001). The basic algorithm we use 
for the radiative transfer is based on the work of Witt (1977). 
The first scattering of the photons is forced to insure that every 
photon contributes to the scattered light (Cashwell & Everett 
1959). The subsequent scatterings are not forced. The direc- 
tion into which the photon is scattered at the nth scattering 
is described by spherical angles, which are determined from 
a Henyey-Greenstein phase function. The scattered light sur- 
face brightness image is constructed using the "peeling-off" 
technique (Yusef-Zadeh et al. 1984). Dust albedo a and asym- 
metry factors g for each emission line are adopted from the 
Milky Way dust model of Draine (2003). 

We assumed a homogeneous dusty medium with \z\ < 100 
pc, |.x | < 500 pc, and \y\ < 500 pc to include the whole H II 
region. The source center is at (x,y,z) = (0,0,0). The observer 



is assumed to be located at a distance of 3 Mpc, appropriate 
to NGC7793 analyzed by Ferguson et al. (1996b), from the 
source center along the z-axis. We considered four dust den- 
sities corresponding to the number densities of hydrogen «h = 
5, 10, 15, and 20 cm" 3 . The Galactic gas-to-dust ratio (dust 
mass per hydrogen nucleon of 1 .87 x 10~ 26 g/H) was assumed 
for calculation of the dust density. The optical depth from the 
source center to the observer is then 0.58, 1.17, 1.76, and 2.35 
at Ha for hydrogen density of 5, 10, 15, and 20 cm -3 , respec- 
tively. The cloud thickness of 100 pc in the z-direction was 
chosen to have an optical depth at Ha of about unity for the 
hydrogen density of 10 cm" 3 . The spatial extents of the dust 
cloud in the x- and y-axes were chosen to cover sufficiently 
large angular sizes of the dust-scattered halos. Note that the 
distance to the source and dust density configuration does not 
alter the main results discussed here. 

We consider two simple types of Ha emissivity profile 
e(x,y,z) to obtain general properties of dust-scattered halos, 
which are independent of the detailed shape of the emissivity 
distribution of the photon source. For the first type, Ha pho- 
tons are uniformly emitted from a sphere with a radius r out , 
e.g., e(x,y,z) = constant in r < r out and e{x,y,z) = otherwise. 
In the second type of model, Ha photons are uniformly emit- 
ted from a spherical shell with an outer radius r out and a width 
of Ar = 5 pc, e.g., e(x, y,z) = constant in r m < r < r out and 
e(x,y,z) = otherwise. Here, r m = r out - Ar is the inner radius 
of the spherical shell. In both types of models, the outer radius 
of the source is varied as 10, 20, and 40 pc. The radial coor- 
dinates of Ha photons for sphere models are randomly deter- 
mined by the formula r = r oal £}l 3 , where £ is a uniform ran- 
dom variable between and 1, to ensure uniform emissivity. 
For the shell models, the radial coordinates of Ha photons are 
determined by the formula r = [(r 3 ut -r 3 n )£+r 3 n ] 1//3 . Azimuthal 
and polar coordinates (tp and 8) of photons are isotropically 
determined. Photons are then isotropically emitted from the 
locations (r,9,cf>). 

2.2. General Properties 

We first describe general properties of the dust-scattered ha- 
los. The surface brightnesses of the Ha source should obvi- 
ously be proportional to the total luminosity of the source. In 
the calculations, we assumed that all sources have the same 
total luminosity while the source size and emissivity distri- 
bution shape vary. Therefore, the surface brightness (or in- 
tensity) of the source decreases with the angular size of the 
source when there is no dust-extinction effect. In general, the 
surface brightness would decrease with ~ 6~~ of , where 8 pTO f 
is the angular size of the source emissivity profile convolved 
with the telescope point spread function (including the effect 
of binning brightness profile). 

When the dust radiative transfer effects are considered, the 
total surface brightness (including direct and scattered lights) 
in the source region would still decrease as 6* pro f increases, al- 
though the effects of dust scattering and absorption may alter 
the functional dependence of the brightness on 8 plo f. Figure 
1, in which profiles are binned with a pixel size of 0.29" (cor- 
responding to 4.2 pc at 3 Mpc), shows the resulting surface 
brightness profiles for various source sizes and shapes after 
the dust radiative transfer effect was taken into account. In 
the figure, hydrogen density of 10 cm" 3 was assumed. Solid 
and dashed lines represent the total (direct + scattered) and 
scattered surface brightnesses, respectively. For other hydro- 
gen densities, we obtained qualitatively similar results as in 
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Figure 1. The "surface brightness profile" I(r) in this paper 
represents a horizontal or vertical cut of the two-dimensional 
surface brightness distribution across the source center. 

Surface brightnesses normalized to the same luminosity 
are shown in Figures 1(a) and (b) for sphere and shell mod- 
els, respectively. We note that the surface brightness of the 
dust-scattered halo surrounding a larger source declines less 
steeply near the source boundary. The most interesting point 
is that at distances far from the sources the brightness profiles 
of the dust-scattered halos for various sources are essentially 
the same. This is because an extended source could be ap- 
proximated by a point source, regardless of its shape and size, 
at locations sufficiently far from the source. Dust grains at 
these locations would scatter off the photons as if the photons 
had originated from a point source. Therefore, the bright- 
ness profile of the dust-scattered halo at distances far from 
the source should be independent of the spatial profile of the 
source and depend only on the total luminosity of the source, 
provided that the dust configuration (total optical depth and 
spatial density distribution of dust) was fixed. 

Figures 1(c) and (d) show surface brightness profiles nor- 
malized to the maximum surface brightness for each source 
to compare the relative significance of the dust-scattered halo 
to the brightness of H II region. For a given luminosity, the 
surface brightness in the source region decreases as the source 
size increases, while the brightness profile of the halo remains 
relatively invariable. Therefore, the relative brightness of the 
scattered halo increases with the source size, as shown in the 
figure. We should also note that the relative contribution of 
dust-scattered components to the total intensity at the source 
region increases with source size. In Figure 1 («h = 10 cm -3 ), 
the contribution of scattered light to total (direct + scattered) 
surface brightness at the center is ~ 13-46%. The fraction of 
the scattered flux integrated over all regions (H II + halo re- 
gions) to the total (direct + scattered) flux over all regions is 
about -~62-63%. This ratio is more or less constant, regard- 
less of the source size and shape, for a given optical depth. 
The ratio of the scattered flux integrated over the halo region 
to the total flux integrated over all regions is found to be about 
^42-57%. The ratio of scattered flux in the halo region to the 
total flux over all regions tends to decrease slightly with the 
source size. This is because the scattered fraction in the H II 
region increases with the source size, while the total scattered- 
fraction is independent of the size. The scattered fractions 
measured in various ways are summarized in Table 1. Ob- 
viously, the contribution of dust-scattered light to total flux 
depends on the dust optical depth as well as the geometrical 
configuration of the dust cloud. The dust-scattered compo- 
nent becomes more important when the dust opacity (here, 
hydrogen density) increases. 

We also note that the surface brightness distribution I(r) of 
the individual dust-scattered halos at a large distance r > ro « 
l m f p can be represented well by the following formula: 

I(r) = I(r )(^) 2 e- (r - r ° )/ '\ (1) 

Here, /(ro) and Z m f p = («HC e xt) _1 are the surface brightness at 
r = ro and the mean free path, respectively. a ext is the dust 
extinction cross-section per hydrogen. Figure 2 shows the 
r 2 /(r) profiles of the dust-scattered halos surrounding a spher- 
ical source with a radius of 10 pc, in linear-logarithmic scale, 
in media with various dust densities corresponding to n H = 2, 
5, 10, and 20 cm -3 , which are typical in the ISM. We fitted the 



asymptotic profiles at r > Z m f p , and obtained the best-fit values 
of the characteristic scale length r c , which are consistent with 
the mean free path within a factor of 1 .4. The mean free paths 
and the best-fit scale lengths are shown in the figure. Higher 
density resulted in a more extended profile than that expected 
by Zmfp due to multiple scattering. The formula will be used 
to understand the increase of the [N II]/Ha and [S II]/Ha line 
ratios with decreasing Ha intensity in Section 4 as well as the 
global Ha morphology in face-on galaxies in Section 5. 

Note that our results were inferred from two different spa- 
tial distributions of the photon source, i.e., sphere and shell, 
and therefore they are independent of the spatial distribution 
of Ha emission in H II regions. We confirmed the results 
again using more realistic emissivity distributions predicted 
from photoionization models in Section 3. We also note that 
scattering by dust within H II regions is not a major concern of 
the present study. The halos, which we are mainly interested 
in, are the results of scattering by dust in the lines of sight 
passing through the ISM outside of H II regions. Dust within 
H II regions contributes to the scattered light that is observed 
in the sightlines toward H II regions. 

2.3. Comparison with Observations 

Ferguson et al. (1996a) found that there is no evidence for 
halos around bright OB associations in the /?-band continuum 
images, which they assume to be equally extended as the Ha 
halos surrounding H II regions if the Ha halos are indeed the 
result of dust scattering. Therefore, they concluded that the 
role of the dust scattering is negligible in the Ha emission 
surrounding bright H II regions. We here note that the Ha 
sources are mainly concentrated in relatively few bright H II 
regions of substantial Ha luminosity. On the other hand, the 
diffuse /?-band continuum background in galaxy images is far 
more uniformly distributed spatially, because the /?-band con- 
tinuum sources do not just consist of OB associations but are 
made up mostly by later-type stars, including late-type giants, 
which have a much more uniform distribution than that of the 
OB associations. Accordingly, the dust-scattered halos sur- 
rounding OB associations in the /?-band images would be hid- 
den among the diffuse background, composed of widely dis- 
tributed stellar sources and the dust-scattered light produced 
by them. 

Typical Ha luminosities of bright H II regions, for example, 
in NGC 7793 analyzed by Ferguson et al. (1996b) are > 10 39 
erg s" 1 . Assuming a typical Ha luminosity of 2 x 10 39 erg 
s" 1 and case B recombination, Lyc luminosity of a typical 
OB association is Q(H°) = 1.4 x 10 51 photons s" 1 , which is 
equivalent to the luminosity produced by 160 stars of 07V 
type (see Table 3 which will be explained in Section 3). Us- 
ing the Kurucz stellar model (Castelli & Kurucz 2003), R- 
band luminosity of the OB association composed of 160 07V 
stars is 7 x 10 24 erg s" 1 Hz" 1 . We verified the estimated lu- 
minosities with an evolutionary synthesis code Starburst99 
(Leitherer et al. 1999, 2010). Assuming a Salpeter initial 
mass function with a total mass of 10 4 M Q inthe mass range of 
6Mq <M < 100M Q , an 1.3Myr-oldOB association produces 
the same Lyc and /?-band luminosities as the above Kurucz 
model. Adopting the total /?-band magnitude of 8.71 mag, 
corresponding to 1.01 Jy (1 Jy = 10~ 23 erg s" 1 cm" 2 Hz" 1 ), 
and the major and minor diameters of 9.3 and 6.3', taken 
from NASA/IPAC Extragalactic Database (NED), the aver- 
age background surface brightness in /?-band is 2.6 x 10 5 Jy 
sr" 1 . The surface brightness of the galaxy shows a roughly 
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exponential decrease with a scalelength of l'.l (~ 1 kpc) 
(Carignan 1985). A similar brightness to the above aver- 
age value is found at the galactocentric radius of ~ 3' from 
the radial profile. Using these values, we calculated the ex- 
pected surface brightness profiles in /?-band for OB associa- 
tions with radii of 10 and 30 pc, as shown in Figures 3(a) and 
(b), respectively. In the calculations, we assumed that pho- 
ton sources are uniformly distributed within the assumed radii 
and the hydrogen density was 10 cm -3 . The calculated im- 
ages were convolved with a two-dimensional gaussian func- 
tion with a full width at half maximum of 1.2" (17 pc at 3 
Mpc) and the brightness profiles were obtained with a bin size 
of 0.44" to match the seeing and pixel size of the observations 
of Ferguson et al. (1996b). Black curves show the resulting 
surface brightness profiles in the presence of the diffuse back- 
ground, while red curves show the cases where there is no 
constant background. It is now clear that the dust-scattered 
halos of OB associations in /?-band images are not observ- 
able. Red curves indicate that the dust-scattered halos can be 
observed only when there is no diffuse background, as in the 
case of Ha emission originating from relatively few H II re- 
gions. 

We also examined the /?-band image of NGC7793 ob- 
served by the SIRTF Nearby Galaxy Survey (SINGS) 
(Kennicutt et al. 2003; Dale et al. 2007), which were obtained 
with the same telescope as used by Ferguson et al. (1996b). 
Figures 3(c) and (d) show examples of the observed /?-band 
brightness profile around OB associations located at the coor- 
dinates (Aa, AS), denoted in the figures, relative to the galac- 
tic center. In Figure 3(c), of which the samples are located at 
a galactocentric radius of ~ 3', the black curve represents a 
typical sample of OB associations with the highest brightness 
contrast compared to their surrounding background, while the 
red curve shows a case with relatively low contrast. Figure 
3(d) shows two more OB associations that are located closer 
to the galactic center. Because of a higher background level, 
the profiles near the galactic center have lower contrast. In 
general, brightness profiles of OB associations are more or 
less similar to one of the four examples in the figures, after 
multiplying a scale factor depending on location. The profiles 
tend to have sharper contrast at larger galactic radius, but not 
always. We note that the models in Figures 3(a) and (b) repro- 
duce the observed profiles of OB associations in Figure 3(c) 
perfectly. Our models are also able to reproduce the bright- 
ness profiles shown in Figure 3(d) by varying the background 
level and the size of OB associations We, therefore, conclude 
that the presence of a continuum background washes out any 
observable trace of the scattered continuum halo of OB asso- 
ciations. 

The next fact we need to note is that dust clouds far 
from H II regions can scatter Ha photons originating else- 
where in galaxies, as demonstrated for high-latitude clouds 
in the Milky Way (Manila et al. 2007; Lehtinen et al. 2010; 
Witt et al. 2010). Filaments, patches, and arcs, which are of- 
ten observed not only in the Milky Way but also in external 
galaxies, do not necessarily indicate in situ ionized gas. Un- 
like the assertion of Hoopes et al. (1996), these distinct fea- 
tures can be explained by dust scattering as well. In the halos 
of edge-on galaxies, very little correspondence between the 
Ha-emitting filaments and "absorbing" dust structures was 
found (Howk & Savage 2000), which seems to indicate the 
physically distinct origins of dust and Ha-emitting material. 
However, we should note that the scattered intensity toward 



a dust cloud can be higher, lower, or equal to its surround- 
ings (as demonstrated in Mattila et al. 2007), depending on 
the incident intensity from behind of the cloud and the cloud 
optical depth. A cloud with an optical depth > 1 — 2 can ap- 
pear dark, while a cloud with a lower optical depth appears 
bright. Indeed, no correlation and anti-correlation between 
the dust-scattered FUV continuum background and dust has 
been found at high opacity regions (e.g., Seon et al. 201 la). 
As noted by Howk & Savage (2000), the dust filamentary fea- 
tures observed in edge-on galaxies are traced in absorption 
against the background starlight, thereby implying relatively 
opaque dust clouds (Ay > 0.8-2.0). Therefore, no correspon- 
dence between the Ha-emitting filaments and absorbing dust 
features does not necessarily indicate the distinct origins of 
the dust and Ha filaments, but can also be a natural conse- 
quence of the combined effect of absorption and scattering at 
relatively high opacity. 

It is also well known from the observations of face- 
on galaxies that the radial and azimuthal intensity distri- 
butions of the diffuse Ha emission are highly correlated 
with bright H II regions over both small and large scales 
(e.g., Walterbos & Braun 1994; Ferguson et al. 1996a,b; 
Zurita et al. 2002). In fact, the strong correlation of the halo 
brightness with that of the H II region is a natural conse- 
quence of dust scattering. The mean fraction of the diffuse Ha 
emission outside of H II regions is known to be roughly con- 
stant (w 0.5), regardless of Hubble type and the star formation 
rate per unit disk area (Ferguson et al. 1996a,b; Hoopes et al. 
1996; Greenawalt et al. 1998; Zurita et al. 2000; Voges 2006; 
Oey et al. 2007). This mean value matches well with the dust- 
scattered fraction estimated in the present study, as shown in 
Table 1. In fact, it is easy to construct cases where the re- 
flected nebular flux surrounding a star is comparable or even 
larger than the source flux (see Fig. 9 in Witt et al. 1982). 

Zurita et al. (2000) found a weak tendency for the diffuse 
Ha fraction to increase with the galactocentric radius. The 
weak trend can be understood qualitatively, if the diffuse Ha 
emission is dominated by dust scattering. Sizes of H II regions 
would generally decrease with the galactocentric radius (e.g., 
Fich & Blitz 1984). As we noted previously (Table 1), the 
dust-scattered fraction in the halo region slightly increases as 
the source size decreases. Therefore, the scattered diffuse Ha 
fraction may slightly increase with the galactocentric radius. 

Here, we considered Ha halos around single H II regions. 
In a realistic case, the dust-scattered halos from adjacent H II 
regions would merge into each other and produce the perva- 
sive Ha background surrounding H II regions. Therefore, we 
may need more detailed radiative transfer models for individ- 
ual face-on galaxies as a whole system to investigate which 
mechanism, the dust-scattering of Ha or in situ photoioniza- 
tion by Lyc leakage, better explains the overall Ha morpholo- 
gies of face-on galaxies. It is worthwhile to note that some 
quantitative analyses on the Ha morphology were already 
preformed by Zurita et al. (2002) and Seon (2009). In Section 
5, we will show that their analyses for the Ha morphologies 
of face-on galaxies likely better support the dust-scattered Ha 
scenario than the Lyc leakage scenario. 

In summary, the previous arguments based on the spatial ex- 
tent or morphology of Ha halos to rule out the dust-scattering 
origin of the diffuse Ha emission are no longer valid. Instead, 
dust scattering explains the general properties of the diffuse 
Ha morphology very well. 

3. OPTICALLINE RATIOS 
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3.1. Basic Assumptions 

We used the well-known one-dimensional photoionization 
code CLOUDY (version c 10.00 of the code last described 
by Ferland et al. 1998) and the fully three-dimensional Monte 
Carlo photoionization code MOCASSIN (MOnte CArlo Sim- 
ulations of Ionized Nebulae; Ercolano et al. 2003, 2005, 
2008) to model spherically-symmetric radiation-bounded H II 
regions surrounding single ionizing stars in a homogeneous 
medium. We also modeled ionized nebulae in a clumpy 
medium by using the MOCASSIN code. 

Extensive calculations of H II region models were per- 
formed by varying elemental abundances, spectral type of the 
central ionizing source, and hydrogen number density. Vari- 
ous sets of elemental abundances have been previously used 
in modeling the WIM. In the present study, we used six sets 
of abundances, which are summarized in Table 2: new so- 
lar system, B star, Orion nebula, ISM, warm neutral medium 
(WNM) values, and WNM values with enhanced S abundance 
(hereafter, WNM2). The abundance of He is assumed to be 
0.1 of the abundance of H in all six cases. The new solar 
system, ISM, and Orion nebula abundances are as defined in 
CLOUDY. The abundances for B star and WNM are defined 
in Sembach et al. (2000). The WNM abundances represent 
those observed in the well-studied warm diffuse clouds to- 
ward the low-halo star /j, Columbae. The WNM abundances 
were also used by Mathis (2000). The WNM2 abundances 
were used by Wood & Mathis (2004) and are essentially the 
same as the WNM abundances except that the S abundance 
is increased by 0.2 dex. We note that the new solar abun- 
dances adopted in this paper are different from the old solar 
abundances used in Sembach et al. (2000). For the model cal- 
culations using CLOUDY, we assumed Orion dust grains for 
Orion abundances and ISM dust grains for the others, as de- 
fined in CLOUDY. In the MOCASSIN models, no grains were 
assumed. The effect of dust grains is mainly to reduce the ion- 
ized region. However, no significant effect on the line ratios is 
found (Mathis 1986b). We ignored elements other than those 
shown in Table 2. The electron temperature and ionization 
structure of H II regions are mainly determined by H, He, N, 
O, Ne, and S, which are major coolants (e.g., Mathis 1985). 
The elements that were not included in the calculations have 
no significant effect on the models of H II regions. 

Spectra of ionizing stars are also needed to model pho- 
toionized nebulae. We considered only dwarf stars with spec- 
tral types from 03V to B1V. For each spectral type, we 
adopted a stellar temperature, radius, and effective gravity 
derived from the evolutionary tracks in Straizys & Kuriliene 
(1981). We calculated the spectral energy distribution and 
the hydrogen ionization luminosity Q(H°) for each spec- 
tral type with a grid of Kurucz models, which are plane- 
parallel line-blanketed model atmospheres in local thermo- 
dynamic equilibrium (LTE) (Castelli & Kurucz 2003). Ta- 
ble 3 summarizes the stellar temperatures, radii R*, H- 
ionizing (Lyc) luminosities Q(H°), He-ionizing luminosities 
(2(He°), and gravities g t for the spectral types used in this 
study. The non-LTE (NLTE) WMBASIC model atmospheres 
(Pauldrach et al. 2001; Sternberg et al. 2003) are also used 
to calculate the one-dimensional H II regions and the results 
are compared with those obtained with the LTE atmospheres. 
Most of the photoionization models in the present study are 
calculated by using Kurucz's LTE atmospheres, unless other- 
wise specified. 

The model calculations using CLOUDY were terminated 



when the fraction of neutral hydrogen (x e dge) reached 0.99, at 
which point the emissivities are very low. We found no sig- 
nificant difference when further calculations were performed 
beyond this value. Note that most of the one-dimensional 
WIM models used nebulae models with x e d ge = 0.95 (Mathis 
1986a) or combinations of models with x e d ge = 0.10 and 
■Xedge = 0.95 to explain line ratios (Domgtirgen & Mathis 1994; 
Sembach et al. 2000; Mathis 2000). For the models calculated 
with MOCASSIN, the sizes of density grid structures were 
decided to be at least -~20% larger than the Stromgren radius 
calculated with CLOUDY to take the effect of absence of dust 
into account and ensure the models cover the whole ionized 
regions. Reflection symmetry in the x-, y-, and z-directions 
was assumed for fast calculations. The number of grids is 
33 x 33 x 33 for the MOCASSIN models. 

We present only the emission lines Ha, [N II] A6583, and 
[S II] A6716, which were best investigated in the Milky Way, 
and He I A5876, which can best constrain the spectral type 
of ionizing stars. The volume emissivities of Ha, He I, [N II], 
and [S II] lines obtained using the CLOUDY and MOCASSIN 
models were used as input emissivity distributions to cal- 
culate the dust-scattered halos of the optical lines surround- 
ing H II regions. For the dust halo models of spherically- 
symmetric H II regions calculated with the CLOUDY code, 
the initial radial coordinates of photons were determined by 
the volume emissivity e of each optical line at a radius r, i.e., 
P(r)dr = er 2 dr/ J er 2 dr. In the dust-scattering models of the 
H II regions calculated with the MOCASSIN code, the ini- 
tial coordinates {x,y,z) of photons were randomly selected 
to follow the volume emissivity e(x,y,z), i.e., P(x,y,z)dV = 
e(x,y,z)dV / J e(x,y,z)dV, and their initial directions were ran- 
domly chosen from an isotropic distribution. 

3.2. One-dimensional Models 

We first modeled the one-dimensional H II regions by vary- 
ing hydrogen number density «h, elemental abundances, and 
spectral type of the central ionizing source. Hydrogen num- 
ber density was varied as 0.1, 1, 10, and 100 cm -3 . Fig- 
ure 4 shows brightness profiles of line intensities (Ha, He I, 
[N II], and [S II]) and intensity ratios (He I/Ha, [N II]/Ha, and 
[S II]/Ha), projected onto the sky, for various spectral types 
of ionizing source. In the figure, the hydrogen density and el- 
emental abundances were assumed to be 10 cm" 3 and WNM, 
respectively. The curves from the outermost to innermost cor- 
respond to later type stars progressively. The outermost curve, 
corresponding to the largest H II region, is obtained from an 
03V star, and the innermost curve corresponds to a Bl V star. 

Emissivities of the low ionization lines [S II] and [N II] in- 
crease suddenly near the ionization boundary of the H II re- 
gion, for most of ionizing stellar types except Bl V, and then 
drop, whereas the Ha and He I recombination lines (RLs) de- 
crease continuously. In other words, both [N II] and [S II] arise 
primarily in the transition region between the ionized and par- 
tially ionized zones of the H II region (e.g., Evans & Dopita 
1985). Therefore, the line ratios of [S II]/Ha and [N II]/Ha 
increase rapidly at the ionization boundary. Such a rapid in- 
crease in the line ratios at the ionization boundary is shown 
in the Orion nebula (Pogge et al. 1992; Sanchez et al. 2007; 
Mesa-Delgado et al. 201 1). The giant H II region NGC 595 in 
M 33 also shows an increase in the line ratios at the ionization 
boundary (Relano et al. 2010). 

The line ratio [N II]/Ha shows a less rapid increase at the 
boundaries than [S II]/Ha. S + is much more confined to the 
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outer regions of the Strtimgren sphere than is N + (Mathis 
2000). The reason is that the change in the predominant 
ionization state between the central and outer H II regions is 
smaller for N because of its higher ionization potential for N + 
— > N 2+ : the ionization potentials for N + — > N 2+ and S + — > S 2+ 
are 29.6 eV and 23.3 eV, respectively, and therefore S exists 
in the doubly ionized state out to larger radius than N. 

For the hottest (03 V and 04V) stars, which emit a sufficient 
number of photons at higher energies than the He ionization 
potential (24.6 eV), photons with E < 24.6 eV run out slightly 
faster than those with E > 24.6 eV at the ionization boundary. 
This resulted in a less rapid decline of the He I line intensity 
than Ha and a weak spike in the line ratio He I/Ha at the 
ionization boundary, as shown in Figure 4. However, for later 
stars, He becomes neutral at the inner zone relative to H, and 
thereby He I/Ha line declines rapidly at a much smaller radius 
than the Stromgren radius. 

As already noted, the surface brightness of the dust- 
scattered halo depends on the total luminosity of the source 
only, for a fixed dust configuration. Therefore, the intensity 
ratio of a pair of emission lines in the dust-scattered halo is 
expected to be determined by the ratio of total luminosities 
of the lines, rather than by the ratio of surface brightnesses at 
the main part of the H II region, if their wavelengths are close 
enough to each other. If their wavelengths are different to the 
extent that dust extinction cross sections and albedos are sig- 
nificantly different, their ratio in the dust-scattered halo would 
also differ from the luminosity ratio. In the brightness pro- 
files, the highest [N II] and [S II] intensities occur in relatively 
small regions at the ionization boundary. However, the tran- 
sition zones can contribute significantly to total luminosities 
because the luminosities are defined by integrals 4tt J er 2 dr 
and the highest emissivities of [N II] and [S II] are found at 
the largest radii. Therefore, the line ratios of [S II]/Ha and 
[N II]/Ha at the dust- scattered halo should be larger than the 
values at the central part of the H II region, but lower than the 
values at the boundary as a result of dust-scattering. 

Figure 5 shows the dust-scattered halo profiles of emission 
lines originating from H II regions surrounding 04V and 09 V 
stars, in the homogeneous dusty medium with hydrogen den- 
sity of 10 cm -3 and the WNM abundances. The top panels of 
the figures show brightness profiles of line intensities of Ha, 
He I, [N II], and [S II]. Total (scattered + direct), direct, and 
scattered intensity profiles are denoted by black, red, and blue 
curves, respectively. By comparing the figures, we can con- 
firm the result that the more extended source (04V) produces 
a brighter dust-scattered halo. The dust-scattered component 
within the H II region surrounding an 04V star is also higher 
than the less-extended H II region surrounding the 09V star. 
The bottom panels show the line ratios He I/Ha, [N II]/Ha, 
and [S II]/Ha. It is clear that the line ratios [N II]/Ha and 
[S II]/Ha outside of H II regions are higher than those at the 
central regions and lower than those at ionization boundaries. 
The He I/Ha profile for the 04V star shows a weak spike at 
the ionization boundary. On the other hand, the profile for the 
09V star shows a large dip at the zone where helium is neutral 
but hydrogen is still ionized. In the zone, the directly-escaped 
lines without interaction with dust grains dominate the line ra- 
tio He I/Ha while the dust-scattered light determines the line 
ratio outside of the H II region. These line ratios found in 
the dust-scattered halos outside of the H II regions match ex- 
actly the luminosity ratios found in the H II regions, except 
for a negligible difference in the He I/Ha ratio caused by their 



wavelength difference. Therefore, we can refer to the lumi- 
nosity ratios, for the purpose of estimating the line ratios in 
the dust-scattered halos, without detailed calculations of the 
dust scattering effect. 

Figure 6 shows the luminosity ratios of He I/Ha, [N II]/Ha, 
[S II]/Ha, and [S II]/[N II] for various hydrogen density, abun- 
dances, and spectral types of ionizing stars. In the figure, 
the results obtained not only with the Kurucz (denoted by 
solid lines), but also with WMBASIC (dashed lines with dot 
symbols) stellar models are shown. The line ratio He I/Ha 
is quite independent of elemental abundances and hydrogen 
density, while [N II]/Ha and [S II]/Ha show large variations. 
Early type stars always produce higher He I/Ha ratios than 
late type stars. The Orion abundances yield slightly lower 
He I/Ha ratios compared to other abundances, especially for 
early type stars. The stellar type giving the highest [N II]/Ha 
ratio varies with hydrogen density, for a given set of abun- 
dances. At high densities, the highest [N II]/Ha ratio is ob- 
tained from relatively late type stars as compared to at low 
densities. For example, for the WNM abundances, the high- 
est value at «h = 100 cm" 3 is obtained from an 09 V star while 
an 05 V star gives the highest value at «h = 0.1 cm" 3 . Ignoring 
the Bl V star, the line ratio [S II]/Ha increases with decreasing 
stellar temperature, except for B stars and solar abundances, 
for a given set of abundances. We also note that the low- 
est value of [S II]/[N II] has similar dependence on the stellar 
type as the highest value of [N II]/Ha do. At high densities, 
the lowest value of [S II]/[N II] is obtained from relatively late 
type stars as compared to at low densities. It is also clear 
that abundances affect the line ratios [N II]/Ha and [S II]/Ha 
significantly. The highest [N II]/Ha and [S II]/Ha ratios are 
obtained with the WNM and ISM abundances, respectively, 
for a given stellar type. Their lowest values are generally ob- 
tained with the solar abundances. The abundances of B stars 
also give relatively low values of [N II]/Ha and [S II]/Ha. It 
is also worthwhile noting that the variations of line ratios due 
to the difference in the adopted stellar models is less sensitive 
than due to the abundances difference. 

We should note that the line ratios are intricately tangled 
with the adopted abundances. For example, the S abundance 
of WNM2 is 0.2 dex larger than that of WNM, while abun- 
dances of other elements are the same. The [S II]/Ha line 
ratio obtained from the WNM2 abundances is then obviously 
higher than that from the WNM abundances. However, the 
[N II]/Ha line ratio of the WNM abundances does not remain 
constant, but becomes lower. Note also that even though the N 
abundance between ISM and Orion differs only by 0.05 dex, 
the [N II]/Ha line ratio obtained from the ISM abundances is 
lower than that from the Orion abundances by a factor of up to 
~ 2 because the S abundance of ISM is 0.5 dex larger than that 
of Orion. This is because increasing the S abundance changes 
the heating and cooling of ionized gas in a complex manner. 

We now compare the results with the observed line ra- 
tios. The typical values of [N II]/Ha and [S II]/Ha toward 
local diffuse gas in the Milky Way are sa 0.3-0.6 and ps 0.2- 
0.4, respectively, and they vary substantially from sightline to 
sightline. In contrast, classical H II regions all cluster near 
[N II]/Ha ps 0.25 and [S II]/Ha ps 0.1. The He I line is much 
fainter (He I/Ha < 0.02) and has been studied only in a few 
selected directions. The observed values of [N II]/Ha and 
[S II]/Ha toward local diffuse gas are best represented by the 
models with the WNM and WNM2 abundances. Although 
early types can reproduce the observed values of [N II]/Ha 
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at low densities, they cannot reproduce the observed values 
of [S II]/Ha and He I/Ha. In particular, the observed ratio of 
He I/Ha strongly constrains the spectral type of the ionizing 
stars to 08 or later. The line ratio of [S II]/[N II] w 0.5-0.7 is 
also reasonably well explained with late O stars and the abun- 
dances of WNM and WNM2. We therefore conclude that late 
OB stars are the best candidates to match the observed line 
ratios of He I/Ha, [N II]/Ha, [S II]/Ha, and [S II]/[N II]. In 
addition, our results indicate that the elemental abundances 
should be close to those of WNM and WNM2 to reproduce 
the observed line ratios. The solar and B star abundances do 
not match well the observed line ratios, except for the case of 
the highest density («h = 100 cm" 3 ). The Orion abundances 
match the observed [N II]/Ha ratio, but yield slightly weaker 
[S II]/Ha ratios than the observed values. Inversely, the ISM 
abundances explain the observed [S II]/Ha, but underpredict 
the [N II]/Ha ratio. 

3.3. Three-dimensional Models 

In the previous section, we described the photoionization 
models in the homogeneous ISM. However, the real ISM has 
a filamentary and/or clumpy structure. A photoionized neb- 
ula in a clumpy medium has different ionic abundances and 
temperature variation from the models in a uniform medium 
(Mathis & Wood 2005). The emission coefficients of colli- 
sionally excited lines (CELs), such as [N II] and [S II], in- 
crease strongly with increasing temperature, while those of 
RLs, produced by electron cascades following recombination, 
are not particularly temperature sensitive. Therefore, the lu- 
minosity ratios between CELs and the RL Ha in clumpy mod- 
els would be different from those in uniform models. 

To investigate the effect of dumpiness, we first calculated 
uniform density models (denoted with a letter "U") with the 
MOCASSIN code. The hydrogen density of 10 cm" 3 was 
assumed for the uniform models. We next calculated H II 
models in a clumpy medium, assuming a two-phase medium 
comprising high density clumps embedded in a low density 
interclump medium. The volume-averaged density of hydro- 
gen was also assumed to be 10 cm" 3 , as in the uniform mod- 
els. Two values of dumpiness, which can be defined by the 
density ratio between dense clumps and interclump medium, 
were considered. In the first type of clumpy models ("CI"), 
dense clumps (cells) have a hydrogen density of 20 cm" 3 and 
occupy a volume fraction of 1/3. The remaining interclump 
cells have a density of 5 cm" 3 . In the second type of clumpy 
models ("C2"), clumps have a density of 30 cm" 3 and occupy 
1/4 of the total volume. The hydrogen density of the inter- 
clump medium is then assumed to be 10/3 cm" 3 . The second 
type is clumpier than the first. For each type of model, clump 
cells with high density were randomly cast to occupy the vol- 
ume fraction and the remaining cells were then assumed to 
have lower density. Spectral types of ionizing source and el- 
emental abundances were varied as in the one-dimensional 
models, for each dumpiness (Tables 2 and 3). We note that 
Mathis & Wood (2005) assumed a hierarchical density struc- 
ture for a clumpy medium and compared the resulting ionic 
abundances, but not line ratios, with those from the model 
with uniform density. 

Figure 7 summarizes the luminosity ratios of He I/Ha, 
[N II]/Ha and [S II]/Ha obtained from the uniform and 
clumpy models. We first compare the uniform models with 
the CLOUDY models. The overall dependences of the line 
ratios in the uniform model on the stellar type are generally 



consistent with the results from the CLOUDY model. We 
note, however, that [N II]/Ha ratios are generally higher than 
those from the corresponding CLOUDY models, except for 
the case of the ISM abundances. For early-type stars, the 
[S II]/Ha ratios are slightly lower than those obtained from 
the CLOUDY models. On the other hand, for the late-type 
stars, the line ratios [S II]/Ha are higher than those obtained 
from the CLOUDY models. 

In the figure, the most important result is that the photoion- 
ized H II regions in the clumpy medium produce higher line 
ratios of [N II]/Ha and [S II]/Ha than in the uniform medium. 
The line ratios [N II]/Ha and [S II]/Ha increase as dumpiness 
increases, while the change in the He I/Ha ratio is negligible. 
This is because dense clumps tend to absorb more H- and He- 
ionizing photons than the interclump medium and thus the 
clumps have higher temperature than the gas in the uniform 
models. Luminosities of [N II] and [S II] are strongly biased 
to high density regions because their emission measures are 
highly sensitive to temperature and density. Therefore, their 
luminosities increase with dumpiness while the Ha and He I 
luminosities remain constant within 2%. Another interesting 
result is that the ratios of [N II]/Ha and [S II]/Ha increase 
rather rapidly with dumpiness for early-type stars, but less 
rapidly for late-type stars. The [N II]/Ha ratio of the "C2" 
models is larger than that of the uniform models by a factor of 
- 1 .9-2.0 for 03V star and a factor of - 1 . 1-1 .3 for 09V star, 
depending on the adopted abundances. The [S II]/Ha ratio 
is increased by a factor of ^2.4-2.7 for 03V star and a fac- 
tor of ~ 1.8-2.0 for 09V star. It is also worth noting that the 
[S II]/[N II] ratio becomes less dependent on the stellar type 
as dumpiness increases. 

Note that the line ratios of [N II]/Ha, [S II]/Ha and 
[S II]/[N II] depend strongly on the elemental abundances. 
The line ratio of [N II]/Ha obtained from the "C2" models, 
except those with the ISM abundances, matches well the ob- 
served ratio in the diffuse ISM, regardless of the stellar type. 
In particular, a [N II]/Ha line ratio of up to ~ 0.7 is predicted 
for the WNM abundances and 09V star. However, the stel- 
lar type of the ionizing source is strongly constrained to 08V 
or later-type stars by the observed ratio of He I/Ha, as in the 
CLOUDY models. 

The differences of the [N II]/Ha and [S II]/Ha line ratios 
in the MOCASSIN models from those of the CLOUDY mod- 
els might be due to atomic data differences adopted in the 
codes. For example, the data set for collision strengths, tran- 
sition probabilities, and energy levels was adopted from ver- 
sion 5.2 of the CHIANTI atomic database (Landi et al. 2006). 
However, the CHIANTI database is not included by default 
in calculations with the CLOUDY code. We also note that 
reliable data for the recombination rate coefficients for S are 
not available, except S 2+ . The fact that the models with the 
WNM and WNM2 abundances produced significantly differ- 
ent [N II]/Ha ratios demonstrates the complicated relationship 
between the atomic species. Using the WNM2 abundances in 
the MOCASSIN models, we also obtained the same trend as 
in the case of the CLOUDY code. Therefore, judging which 
models are more reliable over others without reliable atomic 
data sets would be impractical. However, it is true that clumpy 
models produce higher line ratios than uniform models. The 
clumpy models considered in the present study demonstrate 
the effect of dumpiness. A more realistic clumpy medium, 
such as a hierarchical density structure, might help to under- 
stand the effect of dumpiness in photoionization models. 

We now examine the detailed profiles of the dust-scattered 
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halos surrounding H II regions calculated with the clumpy 
models. For the dust radiative transfer models, we assumed 
a dusty cloud with a constant density of «h = 10 cm -3 . The 
physical dimension of the cloud was the same as in the 
CLOUDY models. For the clumpy medium models, we may 
need to use the same clumpy density distribution as used in the 
photoionization models. However, because the optical depths 
within the H II regions are small and the dumpiness in the dust 
cloud does not affect the line ratios, we ignored the dumpi- 
ness in the dust cloud. Figure 8 shows the brightness profiles 
of the [N II]/Ha and [S II]/Ha line ratios in the H II regions 
produced by 04V and 09V stars and their dust-scattered ha- 
los. The results were obtained from the uniform and clumpy 
"C2" models with the WNM abundances. It is obvious that 
both line ratios are increased in the clumpy models. How- 
ever, the maximum ratios show somewhat different behavior. 
The maximum value of [N II]/Ha near the ionization bound- 
ary increased with dumpiness for early-type stars, such as 
03V, while the maximum value decreased for late-type stars. 
On the other hand, the maximum value of [S II]/Ha always 
increased with dumpiness, regardless of stellar type. 

This property may be largely related to the fate of He- 
ionizing photons. As noted in Mathis & Wood (2005), the 
most important factor to determine the ionization states in a 
clumped density distribution is the fate of He-ionizing radia- 
tion, depending on the hardness of the stellar radiation. Hot 
early-type stars produce such plentiful He-ionizing photons 
that the He ionization zone coincides with the H ionization 
zone in a uniform medium. The effect of the clumpy medium 
is to decrease the ionic fraction (He + /He) and reduce the size 
of the He ionization zone. The ionic fraction (N 2+ /N) fol- 
lows a similar pattern as He + . Therefore, the ionic fraction 
(N + /N) increases near the ionization front and thus increases 
the maximum value of [N II]/Ha. However, the He ionization 
zone for cool later-type stars is smaller than the H ionization 
zone. In this case, the He-ionizing photons tend to penetrate 
dilute interclump regions beyond the He ionization zone that 
is expected from uniform models. The same effect is applied 
to N 2+ , and thus (N + /N) becomes a bit lower at the ioniza- 
tion boundary compared to the uniform models. On the other 
hand, since the ionization potential for S + — > S 2+ is smaller 
than the He-ionizing potential and the ionization potential for 
S° — > S + is smaller than the H-ionizing potential, the total 
ionic fraction of S + and the maximum value of [S II]/Ha al- 
ways increase with dumpiness, regardless of the stellar type 
of the ionization source. 

Figure 9 shows sample images of the line ratios calculated 
for an 09V star to demonstrate the difference between the 
uniform and clumpy models. It is clear that transition regions 
in the line ratio maps near the ionization boundary become 
less sharp for the clumpy models. In some cases, the dumpi- 
ness may produce a smooth transition zone and the line ratios 
would smoothly change to the background values, as in the 
model of the H II region surrounding the 09.5 V star ( Oph 
(Wood et al. 2005). 

4. EFFECT OF STELLAR Ha ABSORPTION LINE 

In the previous section, we reproduced the typical values of 
the [N II]/Ha and [S II]/Ha line ratios in the diffuse ISM re- 
gions by photoionization due to late OB stars. However, the 
highest line ratios may require additional non-ionizing heat- 
ing sources, which will be discussed in Section 5. We here 
consider another possibility to explain the elevated [N II]/Hq 
and [S II]/Ha line ratios in the faint Ha regions, which has 



never been discussed sufficiently in analyzing the diffuse op- 
tical emissions. 

4.1. Implications from Integrated Spectra of Galaxies 

We should note that most, if not all, of the spiral galaxies, 
including the Milky Way, showing the diffuse Ha emission 
outside of bright H II regions are Sb or later types, including 
irregulars (e.g., Rossa & Dettmar 2003; Hunter & Gallagher 
1990; Kennicuttet al. 1995), in which the underlying stellar 
continuum is dominated by B- and A-type stars (Kennicutt 
1992a,b). Balmer absorption lines in the underlying stel- 
lar continuum can give rise to underestimation of the diffuse 
Ha intensity and overestimation of the line ratios of forbid- 
den lines, such as [N II] and [S II], over Ha. The equivalent 
widths (EWs) of the Balmer absorption line as a function of 
stellar spectral type are given by Pickles (1998). The max- 
imum Ha absorption EW of ~ 10A occurs at early A-type 
stars. Bica & Alloin (1986) shows the EW of the Ha absorp- 
tion line as a function of age in an evolving stellar popula- 
tion. For a Sb or Sc disk galaxy with a mixed population, it 
appears that an Ha absorption EW of w 4-5A should be ex- 
pected for the background faint star population outside of OB 
associations. The EWs of nebular emission lines in the inte- 
grated spectra of Sb and Sc galaxies are ~ 5- 100A (Kennicutt 
1992a,b; Sodre & Stasiriska 1999). Therefore, the underlying 
Ha absorption line would appreciably increase the line ratios 
[N II]/Ha and [S II]/Ha not only in the integrated spectra, but 
also in the diffuse ISM regions. For example, for galaxies 
with the emission EW of ~ 10A, the line ratios will be in- 
creased by a factor of ~2. The effect of Ha absorption be- 
comes less important for the galaxies with older stellar popu- 
lation. However, even the oldest stellar population with an age 
of > 13 Gyr has an absorption EW of > 1A (Bica & Alloin 
1986). 

Since about half of the Ha emission is known to origi- 
nate from the diffuse ISM regions (DIG) and the mean ratios 
of [N II]/Ha and [S II]/Ha in the DIG is higher than those 
in H II regions by a factor of w 2, the line ratios in the in- 
tegrated spectrum over an entire galaxy are expected to be 
higher than the mean values of H II regions by a factor of 
w 1.5. We therefore expect that the line ratios of integrated 
spectra of star forming galaxies, after correction of the un- 
derlying stellar absorptions, would be significantly different 
from those of individual H II regions in the diagnostic line- 
ratio diagrams (i.e., the BPT diagram; Baldwin et al. 1981; 
Vielleux & Osterbrock 1987), if the enhancement of the line 
ratios in the DIG were not due to the underlying Balmer ab- 
sorption lines. 

Lehnert & Heckman (1994) used the integrated optical 
spectra of galaxies, published by Kennicutt (1992a,b), with 
Ha emission-line EWs > 30A in order to avoid errors in 
the measured line ratios caused by underlying stellar absorp- 
tion lines. They corrected the measured line ratios assuming 
that both the Ha and H/3 stellar absorption lines have EWs 
of 5A and found that in general the line ratios in the inte- 
grated galaxy spectra occupy a region that is intermediate be- 
tween H II regions and the AGNs. The resulting [S II] AA6716, 
6731/Ha ratios, for a given [O III] A5007/H^, were typically 
a factor of ~ 1.5 larger, while the [N II] A6783/Ha ratios were 
more similar to those of H II regions, only being typically 
enhanced by ~ 25%, which is not very significant and well 
within the scatter expected for individual H II regions. They 
therefore concluded that at least 25% of the observed Ha flux 
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in their sample galaxies would arise in the DIG. 

However, more recent results seem to indicate no significant 
difference between individual H II regions and star-forming 
galaxies in the diagnostic diagrams, after careful subtraction 
of the underlying Balmer absorption lines. Using the same 
dataset as Lehnert & Heckman (1994), Sodre & Stasiriska 
(1999) made a better correction for the stellar absorption EWs 
by correlating the observed Ha and H(3 EWs. They found 
that the diagnostic diagram of integrated ratios [O III]/H/? ver- 
sus [N II]/Ha follows the sequence of the giant H II regions 
very well. In their diagram of [O III]/H/3 versus [S II]/Ha, 
the absorption-corrected ratios are also well within the H II 
region sequence (except that those with [O III]/H/3 > 1 lie 
at the upper boundary defined by individual H II regions). 
Moustakas et al. (2006) used a spectral synthesis fitting code 
to find the stellar model continuum, including Balmer absorp- 
tion lines, that optimally reproduces the integrated continuum 
spectrum and estimated pure nebular emission line EWs of 
412 star-forming galaxies after subtracting the best- fit contin- 
uum spectra. They noted that the sequences formed by indi- 
vidual H II regions and star-forming galaxies overlap across 
the full range of emission-line ratios in the [O III] A5007/H/3 
versus [N II]/Ha plane. However, in the [N II]/Ha versus 
[S II]/Ha plane (Fig. 7 of Moustakas & Kennicutt 2006), the 
absorption-corrected line ratios preferentially lie at or slightly 
outside the envelope defined by H II regions. We note that 
this cannot be caused by a relative enhancement or reduction 
of the integrated [N II]/Ha ratios compared to the ratios of H II 
regions, but can be caused by an enhancement of the [S II]/Ha 
ratios in integrated spectra by about 0.1 dex. We also note 
that Fig. 1 of Brinchmann et al. (2004), which was obtained 
by fitting the stellar synthesis model continuum to the Sloan 
Digital Sky Survey data, shows that star-forming galaxies lie 
at the same locations as H II regions in the diagnostic diagram 
of [O III]/H/3 versus [N II]/Ha. 

Unlike the absorption-corrected line ratios in the integrated 
galaxy spectra, the line ratios observed in the DIG are signifi- 
cantly different from individual H II regions in the diagnostic 
diagrams (see Fig. 7 of Moustakas & Kennicutt (2006) and 
Fig. 2(a) of Flores-Fajardo et al. (2009)). Therefore, the en- 
hancement of [N II]/Ha in the diffuse ISM regions may be due 
mainly to the underlying absorption lines. In other words, the 
ratio of [N II]/Ha in the diffuse ISM or DIG regions is not 
likely to be significantly different from that of H II regions if 
the underlying stellar absorption lines are properly appreci- 
ated. In the case of [S II]/Ha, the line ratios observed outside 
H II regions, after correction of the stellar absorption line, may 
be slightly larger than that of H II regions. However, the en- 
hancement would be only by a factor of w 1 .5, taking into ac- 
count the facts that half of the optical lines originate from the 
diffuse regions and the ratio in integrated spectra is only 0. 1 
dex higher than that of H II regions. It is also worth noting that 
the [S II] line arises not only in the H + /H° interface regions of 
ionized nebula, but also in the H°-C + regions of photodis- 
sociation regions (PDRs) illuminated by strong FUV radia- 
tion, because the ionization potential of S° is low (10.36 eV) 
(Petuchowski & Bennett 1995; Storzer & Hollenbach 2000). 
Storzer & Hollenbach (2000) found that the relative strength 
of PDR emission to H + /H° interface emission is « 0.14- 1.5. 
Therefore, in a certain condition with high density and strong 
FUV radiation, the PDR may contribute to the [S II] flux in 
integrated spectra of star forming galaxies. Taking the con- 
tribution of PDR emission into account, the difference in the 



[S II]/Ha line ratio between the DIG and H II regions would 
be reduced even further. However, it is not certain, at this 
point, how significantly the PDR emission would contribute 
to the integrated spectra of star forming galaxies. The [S II] 
flux from PDRs may yield a spatially dependent contribution 
to the [S II]/Ha line ratio because molecular clouds and PDRs 
are rarer, especially, in interarm regions. 

We also note that some authors assumed the same EW for 
all Balmer absorption lines (Ha, H/3, H7, and H<5) (e.g., 
Lehnert & Heckman 1994; Leeetal. 2004). We obtained 
EW(H/3) = 1.52 EW(Ha), EW(H 7 ) = 1.38 EW(Ha), and 
EW(H<5) = 1.36 EW(Ha) by fitting the stellar spectral data 
in Pickles (1998). Therefore, the EWs of H^, H7, and HS 
could be assumed to be the same. However, assuming the 
same EW for Ha and H/3 would cause not only an error in 
dust extinction correction, but also a shift in the diagnostic 
diagrams, such as [O III] A5007/H/3 versus [N II]/Ha. 

In summary, it is likely that the line ratios observed in the 
DIG are not significantly different from those of H II regions. 
In the following section, we will further demonstrate how the 
stellar absorption lines are able to increase substantially the 
line ratios in the diffuse ISM regions. 



4.2. Effect on the Lines Ratios in the Diffuse ISM 

If significant variations occur in stellar populations verti- 
cally as well as across spiral arms, the faint Ha regions would 
be significantly affected by the absorption feature in the un- 
derlying continuum. OB stars are usually concentrated at 
star forming regions in spiral arms, while A- and later-type 
stars are relatively widespread. The scale-height of late-type 
stars is also higher than those of earlier types. Therefore, not 
only the direct starlight, but also the scattered light from late- 
type stars would be much more extended than those from OB 
stars. In our Galaxy, for example, the sky can be divided 
into places where OB stars or A-type stars are dominant, as 
clearly shown in Figure 15 of Seon et al. (201 la). Within and 
near H II regions, the stellar continuum would be dominated 
by OB stars and the Ha emission line would overwhelm the 
underlying Ha absorption feature in the continuum compo- 
nent from later-type stars. At locations far from OB stars, 
where later-type stars dominate the continuum, the Ha emis- 
sion line will rapidly decrease while the underlying Ha ab- 
sorption line is far more uniform. Therefore, the absorption 
feature would carve out the emission line. The effect of the 
Ha absorption feature in the continuum would gradually in- 
crease with the distance from bright H II regions and with a 
decrease of the Ha intensity. Therefore, the stellar absorp- 
tion feature should significantly influence faint areas where 
the continuum is dominated by A- or later-type stars. 

We can quantitatively demonstrate how the underlying Ha 
absorption line can significantly affect the measured line ra- 
tios, if the diffuse optical emission lines are mostly caused 
by dust scattering of photons originating from bright H II re- 
gions. First of all, we assume that the surface brightness of 
the dust-scattered continuum from an OB association /°„ t (r) 
decreases as in Equation 1, while the diffuse continuum back- 
ground 4om surrounding the OB association is dominated by 
late-type stars and thus spatially uniform. The observed line 
ratio [N II]/Ha (or [S II]/Ha) at a distance r (from the OB 
association) larger than the mean free path (ro « Z m f p ) in the 
dust-scattered halo of the OB association can be related to the 
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real emission line ratio (I^n] /'h") by 
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where W^ a and W^™ denote the Ha absorption EW measured 
relative to the diffuse continuum background (^ont) an d me 
emission EW relative to the continuum of the OB associ- 
ation (/°nt( r ))' respectively. In the dust-scattered halo, the 
EWs are constant. The emission line ratio /[n"i]// h ™ W1 U be 
constant in the dust-scattered halo, while the line intensities 
decrease together with the continuum intensity from the OB 
association as distance increases. In general, W^™ of bright 
H II regions in disc galaxies ranges from 100A to 1500A, 
with a mean value of 400A (Bresolin & Kennicutt 1997). As 
noted earlier, the most probable absorption EW would be 
W Ha ~ 4A ( Bica & Alloin 1986). We consider, however, not 
only the most probable W^ a = 4A but also the worst case 
of Wjfa = 1A. From Figures 2 and 3, we can assume, for 
« H = 10 cm" 3 , that ftt/O^o) » 4 at r » l mfp = 57 pc and 
ro/r c = 56.9/72.4= 0.79. Figure 10 shows the resulting en- 
hancement factor of the line ratio as a function of distance for 
various combinations of (W^ a , W^™). As shown in the figure, 
the line ratio increases significantly within a few hundred pc. 
The enhancement factor increases at closer distance when the 
absorption EW is higher and/or the emission EW is lower. 

Note that the forbidden-to-Ha line ratios, especially 
[S II]/Ha, are observed to be enhanced by a factor of up to 8 
(or more) compared to the ratios in H II regions (e.g., Fig. 12 
in Madsen et al. 2006). However, such a large increase by a 
factor of > 3 seems to be rare. In this regard, it is important to 
note that the line ratio increases relatively slowly up to a fac- 
tor of ~ 3 and then rapidly beyond the point. Therefore, the 
enhancement factor > 3 would be found only in a small out- 
ermost part of the halo and be rarely observed. In addition, 
we assumed only a single H II region and its scattered halo. 
In reality, there should be other H II regions. The emission- 
lines (Ha, [N II] and [S II]) originating from the neighboring 
H II regions will be scattered to the halo under consideration. 
The scattered Ha from other H II regions may soften the Ha 
absorption effect and the rapid increase of line ratios at the 
outermost part of the halo. 

In the Wisconsin Ha Mapper (WHAM) survey of the 
Milky Way, the bright stars with my < 5.5 were nagged 
to denote possible contamination due to strong Ha absorp- 
tion features (Haffner et al. 2003). However, there still re- 
mains dust-scattered starlight called DGL (Toller et al. 1987; 
Brandt & Draine 2012) and light from stars with my > 5.5. 
The average EW of the diffuse Ha background of our Galaxy 
is estimated to be 1 1 A (Brandt & Draine 2012), assuming the 
interstellar radiation field of Mathis et al. (1983). In the ex- 
ternal galaxies, the Ha absorption feature in the underlying 
continuum both from unresolved stars and the diffuse dust- 
scattered starlight would hamper the estimation of the line ra- 
tios. To the best of our knowledge, Walterbos & Braun (1994) 
was the only one in which the effect of underlying Balmer ab- 
sorption lines was discussed in investigating line ratios in the 
diffuse ISM. For M31, the integrated Ha flux corresponds 
to an EW of only 4-6 A (Walterbos & Braun 1994). There- 



fore, the underlying Ha absorption line is expected to strongly 
affect the observed line ratios in the diffuse ISM regions. 
Broadband surface photometry of M 3 1 shows typical colors 
of B-V = 0.6-0.8, and U-B = 0.0-0.25, which are differ- 
ent from a pure A-star spectrum. This led Walterbos & Braun 
(1994) to claim that the effect due to the Ha absorption line 
are not likely to be important. They also performed a simple 
test, to reduce the effect of underlying Ha absorption line, by 
subtracting a 10% lower continuum from the Ha filter image 
and found no evidence on the effect in the Ha map. However, 
as noted by them, a more detailed test requires high-resolution 
multicolor imaging, combined with spectroscopy to constrain 
the stellar population and its small scale variations. 

5. DISCUSSION 

In this section, we present additional support for our results 
and argue the necessity of reconsidering the dust-scattering 
origin of the diffuse Ha emission more seriously. In fact, we 
should note that there is no direct evidence, given the results 
presented in the previous sections, supporting the contention 
that most of the diffuse Ha emission originates from in situ 
ionized gas. The previous studies of Mattila et al. (2007), 
Witt et al. (2010), and Seon et al. (2011b) have mainly dealt 
with the dust-scattered Ha photons at high Galactic latitudes. 
The present results strongly suggest the significance of dust 
scattering in the diffuse Ha emission not only at high-latitude 
clouds but also near H II regions. 

5.1. Dust-Scattered Ha Halos 

There is observational evidence that at least some parts of 
the diffuse Ha emission outside of classical bright H II regions 
may be in fact caused by dust scattering. If the dust scatter- 
ing is the dominant source of the diffuse Ha halos around 
the bright H II regions, the halos would not be observable or 
would be very faint in a radio recombination continuum that 
is almost free from dust-scattering. In this way, at least some 
part of the extended Ha halo of the Orion H II region was at- 
tributed to the dust-scattered halo (Subrahmanyan et al. 2001; 
O'Dell & Goss 2009). More direct evidence of dust-scattered 
origin of the Ha halo can be provided by observations of the 
Ha polarization. Topasna (1999) attempted to detect Ha po- 
larizations by selective extinction in the Monoceros supernova 
remnant, the Rosette Nebula, and the North America Nebula. 
While these were not detected, Ha polarizations by scatter- 
ing in dust shells around the Rosette Nebula and the North 
America Nebula were instead observed. 

Additional support for the importance of dust scattering in 
the diffuse Ha emission is provided from the large scale sur- 
veys of our Galaxy. Kutyrev et al. (2001, 2004) carried out 
a pilot survey of the Galactic plane with 2.17 /im Br7 hy- 
drogen RL, which is relatively free from dust-scattering, and 
found a typical filling factor of a Br7-emitting gas of ~ 1 %, 
which is much smaller than the value derived for Ha. Re- 
cent observations using the Wilkinson Microwave Anisotropy 
Probe have found that the ratio of the free-free radio contin- 
uum to Ha is surprisingly low in the WIM (Davies et al. 2006; 
Dobler & Finkbeiner 2008a,b; Dobler et al. 2009; Gold et al. 
201 1). These discrepancies likely suggest the significant role 
of dust scattering in producing the diffuse Ha emission out- 
side of bright H II regions. However, we note that this dis- 
crepancy is not just present in regions where dust scattering 
may contribute, but even in regions that appear to be rather 
dust free, e.g., the Gum nebula. The cause of this observed 
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fact may therefore be more complex and involve other factors 
besides scattering. 

It is also of interest to note that a correlation between Ha 
and the anomalous dust emission at ~ 30 GHz, which is at- 
tributed to spinning very small dust grains, is found in the 
WMAP data (Dobler & Finkbeiner 2008a,b). The spectrum 
of Ha-correlated microwave emission does not follow the 
expected free-free spectrum from the WIM, but the electric 
dipole emission from rapidly rotating dust grains. Since then, 
there have been a number of studies showing the association 
of peaks in the anomalous dust emission with individual in- 
terstellar clouds (e.g., Vidal et al. 2011). The spinning dust 
microwave emission is simply the signal from the smallest 
component of the grain size distribution in interstellar space. 
The correlation of the anomalous dust emission with Ha is 
explained most readily, if a significant fraction of the Ha is 
also coming from dust via scattering. 

The photoionization models of the diffuse Ha emis- 
sion assume no significant absorption of the Lyc 
(Domgorgen & Mathis 1994). A widely accepted view 
for the almost-free escape of Lyc from the bright H II regions 
is to assume a highly clumpy or turbulent density structure 
of ISM with a lot of almost-empty space. Ferguson et al. 
(1996a,b) found a clear correlation between the diffuse Ha 
intensity and the surface brightness of H II regions on both 
large and small scales. This correlation was regarded to 
be evidence that the Lyc leakage is the dominant source 
responsible for producing the diffuse Ha photons. However, 
as we demonstrated, the radiative transfer models for the 
dust-scattered halos around single H II regions can reproduce 
the basic properties of the Ha halos surrounding H II regions 
in face-on galaxies. An advantage of the dust-scattering 
scenario investigated in the present study over the in situ 
ionization scenario would be the fact that the diffuse Ha 
fraction of « 50% can be explained naturally. 

Since the analysis was based on the halos surrounding sin- 
gle H II regions, it is necessary to investigate if dust scattering 
can indeed explain the global Ha morphology, which may be 
an overlapping of the dust-scattered halos surrounding indi- 
vidual H II regions, of the face-on galaxies. In this regard, 
we point out that the analyses of Zurita et al. (2002) and Seon 
(2009), which were initially intended to explain the global Ha 
morphology in terms of photoionization models, can be inter- 
preted as accounting for the morphology by the overlapped 
dust-scattered halos. Zurita et al. (2002) attempted to model 
the diffuse Ha morphology seen in face-on galaxies assum- 
ing that the Lyc leakage from bright H II regions would be 
the predominant source of the diffuse Ha. They assumed that 
a constant fraction of the Lyc leaks from bright H II regions 
and transferred through the two dimensional ISM disk with no 
significant absorption, and claimed that the diffuse Ha sur- 
face brightness of NGC 157 can be reproduced reasonably 
well with this scenario. Seon (2009) extended the analysis 
to three dimensions assuming an exponential disk, and per- 
formed a similar analysis for a face-on galaxy, M 51. We 
then concluded that the effective hydrogen density required 
to explain the diffuse Ha emission with the "standard" sce- 
nario was too small (~ 10~ 5 of the generally known value) 
to be reconciled with the effective density inferred through 
the turbulence properties of ISM. In other words, unlikely ex- 
treme topology of the ISM is required to explain the diffuse 
Ha morphology by photoionization in a turbulent or clumpy 
medium. 



These analyses, however, did not entail self-consistent pho- 
toionization modeling, as noted in Wood et al. (2010). The 
photoionization model should be iteratively performed until 
the ionization fraction and temperature at each points con- 
verge. Instead, the penetration of photons through ISM was 
characterized using an exp(-r) factor, after being diluted by a 
geometric factor 1 / r 2 , where r is the effective optical depth, 
and the photons interacting at a point were assumed to be re- 
emitted isotropically. This radiative transfer process in fact 
describes a single dust-scattering with an arbitrary dust ex- 
tinction cross-section cr e ff and an isotropic scattering phase 
function (asymmetry phase factor g = 0). In our simula- 
tions, it was found that the Ha luminosity from the diffuse 
region (DIG) is approximately half of the total Ha luminos- 
ity originating from H II regions. In the context of a sin- 
gle dust-scattering, this is equivalent to assuming the dust 
albedo of a ~ 0.5. From the analysis, it was also found that 
the effective cross-section per hydrogen nucleus is given by 
(Jeff ~ (2-4) x 10~ 22 cm 2 for a reasonable hydrogen density. 
Note that the obtained effective cross-section is a factor of 
?» 10~ 5 smaller than the photoionization cross-section of Lyc, 
but accords well with the dust-extinction cross-section of Ha 
photons of cr ext = 3.8 x 10" 22 cm 2 /H (Draine 2003). The ab- 
sorption coefficient obtained by Zurita et al. (2002) is also 
consistent with the value. Therefore, the results indicate that 
the diffuse Ha morphologies in the face-on galaxies accord 
well with the dust-scattered halos of Ha photons, contrary to 
the claim of Zurita et al. (2002) that the analysis supports the 
photoionization by Lyc leaked out of bright H II regions. 

In Section 2, we found that the surface brightness of the 
individual dust-scattered halos can be approximated well by 
the formula I(r) oc exp(-r)/r 2 , which is in fact the very ex- 
pression used in Zurita et al. (2002) and Seon (2009). These 
results seem to strongly suggest the dust-scattering origin of 
the diffuse Ha emission in the face-on galaxies. Obviously, 
the radial profile is not consistent with those of the radiation- 
bounded H II regions. However, it is not clear at this stage 
whether or not the radial profile of the photoionized Ha emis- 
sion line that is predicted from the density-bounded models 
accords with I(r) oc exp(-r)/r 2 . We plan to examine the ra- 
dial profiles of density-bounded H II regions in the future. We 
also need to develop more realistic models of the global Ha 
morphologies of face-on galaxies in the contexts of not only 
dust scattering, but also photoionization. 

5.2. Optical Line Ratios 

The [N II]/Ha and [S II]/Ha line ratios are observed to be 
higher in the diffuse Ha regions than in bright H II regions 
(Haffner et al. 2009). This was believed to be the strongest 
evidence against the dust-scattering origin of the diffuse Ha 
emission, because the dust scattering would keep the line ra- 
tios constant. The enhanced optical line ratios outside of 
H II regions have been generally explained by a dilute radi- 
ation field transferred from O stars within H II regions to the 
diffuse ISM (Mathis 1986a; Sokolowski & Bland-Hawthorn 
1991; Domgorgen & Mathis 1994). 

However, we note that the argument to reject dust scattering 
has been based solely on comparisons with the line ratios from 
bright H II regions. If the line ratios from the ionized nebu- 
lae due to the late O- or early B-type stars are compared with 
the typical line ratios in the diffuse Ha regions, the present 
scenario seems to explain the higher line ratios than those in 
bright H II regions equally as well as the previous photoion- 
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ization models do. It was also found that the [S II]/Ha ra- 
tios observed in H II regions around late-type stars tend to 
be higher than those found around early-type stars (Figure 3 
of Reynolds 1988). O'Dell et al. (2011) pointed out a strik- 
ing similarity of the line ratios between Barnard's Loop, the 
Orion-Eridanus Bubble, and the typical WIM samples, which 
can be explained by the photoionized nebulae due to stars with 
temperature of < 35, 000 K corresponding to late O-types (see 
also, Seonetal. 2011b). This also supports the present con- 
clusion. 

Weak [O III] A5007/Ha and He I/Ha emission line ratios in 
the galactic disks also indicate that the spectrum of the diffuse 
interstellar radiation field is significantly softer than that from 
the average Galactic O star population (Reynolds & Tufte 
1995; Madsen et al. 2006). These line ratios provide the 
strongest constraint on the spectral type of the original ion- 
izing source of the diffuse Ha and other optical emission 
lines. However, [O III] A5007/H/3 or [O III] A5007/Ha in 
several halos of edge-on galaxies are found to increase with 
height above the galactic plane (Tullmann & Dettmar 2000; 
Miller & Veilleux 2003). The rise may indicate that other 
mechanisms are needed. 

Seon et al. (201 lb) estimated how a large portion of the ion- 
izing power required for the Ha background in our Galaxy 
can be provided by late OB stars, and concluded that 09 and 
later-type stars can account for at least one-half of the required 
recombination rate of rQ « 2.5 x 10 6 s" 1 cm" 2 . We note that 
the earliest stellar type consistent with the observed He I/Ha 
line ratio is 08 from our analyses, and the fraction of Lyc 
luminosities of 08 and 09 stars in the total Lyc luminosity 
due to O stars is ~ 18% (Table 1 in Terzian 1974). Including 
08 stars, late OB stars can provide the recombination rate of 
(2.3-6.2) x 10 6 s" 1 cm" 2 . It is therefore clear that most of 
the required ionizing power in the Galaxy would be explained 
with late OB stars. In fact, there is a relative lack of early O 
stars in the solar neighborhood, as noted in Brandt & Draine 
(2012). Therefore, the diffuse optical emission lines in our 
Galaxy are consistent with those produced by late OB stars 
and their scattered light from the local ISM. 

In our models using the MOCASSIN code, the typical val- 
ues of the [N II]/Ha and [S II]/Ha ratios were well repro- 
duced with most of the abundances except the ISM abun- 
dances. With the CLOUDY code, abundances close to those 
of WNM were required. Sembach et al. (2000) and Mathis 
(2000) also reproduced the observed line ratios assuming the 
WNM abundances. In our previous paper (Seon et al. 201 lb), 
we mistakenly stated that the elemental abundances required 
to explain the line ratios in the WIM are close to the B star 
abundances both in Sembach et al. (2000) and in O'Dell et al. 
(201 1). This should be read as the WNM abundances instead 
of the B star abundances. 

The NLTE stellar models may produce significantly differ- 
ent amount of Lyc photons from the LTE models. Neverthe- 
less, we obtained consistent results from two different model 
atmospheres (Kurucz or WMBASIC) in Section 3.2. Note 
that it is the combination of the line blanketing and spheri- 
cal geometry, not exclusively the NTLE aspect, that is cru- 
cial in determining the spectral energy distribution (SED) and 
Lyc luminosity of stellar models (Aufdenberg et al. 1998). 
Aufdenberg et al. (1999) also found that a significant differ- 
ence between the spherical and plane-parallel models is found 
only at surface gravities below logg* = 3.5, corresponding 
to giant stars. We also note that the calibration scale of 



stellar parameters (temperature and gravity) as a function of 
spectral type varies when different model atmosphere codes 
are used (Martins et al. 2005; Simon-Diaz & Stasihska 2008). 
Simon-Diaz & Stasihska (2008) investigated the impact of the 
modern model atmospheres (including plane-parallel stellar 
models) of dwarf stars on H II regions. They concluded that 
the predicted SEDs in the energy range of 13-30 eV, which 
are most important in producing Ha, [N II], [S II], and He I 
lines, are in good agreement between the codes when the 
models are compared using the calibration scales appropriate 
to each code. Therefore, our main results obtained by using 
dwarf stars are not sensitive to the adopted atmosphere mod- 
els. Since the number of giant stars is relatively smaller than 
that of dwarfs, it is also unlikely that the main conclusions 
will be significantly altered, even when giants are taken into 
account. Fully three-dimensional simulations with a set of 
NLTE input spectra, including the effect of giant stars, may 
be needed for a more detailed study. 

The results of Topasna (1999) are worth noting in attempt- 
ing to understand the line ratios of the diffuse Ha regions. The 
[S II]/Ha ratio observed in the Rosette Nebula increased from 
< 0.2 to ~ 0.5 with radius out to the boundary of the H II 
region at r < 35', which was identified by polarization and 
4850 MHz radio continuum emission. Beyond the H II region 
boundary (r « 35' -60') where dust scattering is dominant, 
there were no large variations in the ratio as a dust-scattered 
halo would preserve the line ratio at the H II region bound- 
ary (see Fig. 4.27 of Topasna 1999). The ( Oph H II region 
also shows the same trend in the [N II]/Ha and [S II]/Ha ra- 
tios (Fig. 2 of Wood et al. 2005). This trend is exactly what 
is expected from the dust-scattered halo. However, farther out 
beyond these regions with constant ratios, the line ratios in- 
creased to higher values. These higher ratios may indicate the 
existence of other sources close to the nebulae to elevate the 
line ratios, as will be discussed in the following paragraph. 

It is thought that the line ratios increase as the Ha intensity 
decreases. This trend does not appear to be consistent with 
the dust scattering scenario, because dust scattering should 
preserve the ratios. We therefore need to explain the trend 
in the context of the dust scattering scenario. Note that early 
O stars in association are more concentrated at the galactic 
plane. Near an early O star, the dust-scattered Ha halo would 
be dominated by the H II region due to the O star and the Ha 
intensity will be very high. The [N II]/Ha and [S II]/Ha ra- 
tios would be then relatively low in the region. At distances 
far from the O star, the Ha intensity will be dominated by 
the dust-scattered light originating from the H II regions sur- 
rounding late OB stars, which are more abundant than early 
O stars, and the line ratios will then be increased. As we 
demonstrated using the clumpy models, some part of the en- 
hancement or fluctuation in the line ratios may be attributed to 
the variation of dumpiness from sightline to sightline. Note 
that the [S II]/Ha ratio of ~ 0.23 in the C Oph H II region 
(Wood et al. 2005) is close the ratio expected from the CI 
model, while the [S II]/Ha ratio of ~ 0.5 in the Rosette Neb- 
ula (Topasna 1999) needs a clumpier model than the C2 case. 
In LDN 1780, the ratio of ~ 0.16 (Witt et al. 2010) is some- 
where in the range between the U and CI model. However, 
the highest line ratios may need other explanations. The high- 
est line ratios can result from non-ionizing heating sources, as 
discussed in Seon et al. (2011b), and/or from the underlying 
Balmer absorption line in the stellar continuum of relatively 
late-type stars, as discussed in Section 4. 

The heating sources to elevate the line ratios may in- 
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elude shocks, photoelectric heating, turbulent mixing 
layers, and/or galactic fountain gas Reynolds & Cox 
(1992); Slavinetal. (1993); Collins & Rand (2001); 
Raymond (1992); Shapiro & Benjamin (1993). Recently, 
Flores-Fajardo et al. (2011) proposed hot low-mass evolved 
stars (HOLMES), which may be plentiful in the thick disks 
and lower halos of galaxies, to explain the observed increase 
of [O m]/H/3, [O II]/H/3 and [N II]/Ha with increasing 
distance to the galactic plane in an edge-on galaxy, NGC 891. 
We also note that the H II region associated with the 
B0.5+sdO binary Per is elevated in [S II]/Ha (Madsen et al. 
2006). These stars produce a much harder radiation field 
than massive OB stars, and are therefore able to produce 
H II regions with high electron temperatures. It is possible 
that all these supplementary sources may play some roles 
in enhancing the line ratios and their contributions might 
significantly vary from sightline to sightline. The optical 
emission lines due to these sources providing the highest line 
ratios should also be scattered into more extended regions 
than those originally produced by the sources. It should be 
noted that the conventional photoionization models by O 
stars in the galactic plane also need these heating sources to 
explain the highest ratios (Reynolds et al. 1999). Therefore, 
the present dust scattering scenario can be an alternative 
explanation for the diffuse optical Ha emission that is equally 
as plausible as the in situ ionization scenario. 

In Section 4, we argued that the [N II]/Ha line ratio (prob- 
ably, [S II]/Ha as well) in the diffuse ISM regions seems to 
be consistent with that in bright H II regions, if the underlying 
Balmer absorption lines are taken into account. It was also 
demonstrated that the line ratios should indeed increase at the 
faint Ha regions outside bright H II regions, if the diffuse Ha 
emission is predominantly caused by dust scattering. We note 
that a similar argument to explain the anti-correlation of the 
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line ratios with Ha intensity can be applied even to the in 
situ photoionization scenario. Then, the non-ionizing heating 
sources proposed to explain the highest line ratios may not 
be required anymore, regardless of what would be the main 
origin of the diffuse Ha emission, dust-scattered light or in- 
situ ionized gas. However, the dust scattering origin is more 
preferable if this is indeed the case. 



6. SUMMARY 

We demonstrated that the basic morphological properties 
of the diffuse Ha emission can be well explained by the dust- 
scattered halo surrounding the H II nebula, as opposed to the 
previous belief that dust scattering does not accord with the 
Ha and /?-band observations. We also found that the optical 
line ratios of He I/Ha, [N II]/Ha and [S II]/Ha observed in the 
diffuse ISM outside of bright H II regions can be well repro- 
duced by the dust-scattered halo scenario, wherein the optical 
lines originate from H II regions ionized by late O- or early 
B-type stars in the media with abundances close to WNM and 
are scattered off by the interstellar dust. It was also shown 
that the predicted [N II]/Ha and [S II]/Ha line ratios increase 
with the dumpiness of ISM. We also showed that the Balmer 
absorption lines in the underlying stellar continuum seem to 
explain the rise of line ratios in the faint Ha regions. We 
therefore conclude that the dust-scattering origin of the dif- 
fuse optical emission lines cannot be ruled out. Instead, the 
dust scattering can explain the observed properties of the dif- 
fuse optical emission lines very well. 
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tic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract 
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TABLE 1 

Scattered fractions for various models 



model 


scattered intensity / total intensity 
at source center 




scattered flux of halo 
/ total flux over all regions 


scattered flux over all regions 
/ total flux over all regions 




5 cm-' 


10 cnr J 


15 cnr J 


20 cnT J 


5 cm 1 


10 cm -3 


15 cnT 3 


20 cm -3 


5 cm 3 


10 cm -3 


15 cm' 3 


20 car 3 


sphere 10 pc 


6.93% 


12.5% 


18.4% 


24.1% 


41.2% 


56.9% 


65.5% 


71.2% 


45.5% 


63.0% 


72.7% 


79.0% 


sphere 20 pc 


11.5% 


21.9% 


31.1% 


39.4% 


37.5% 


51.1% 


58.3% 


62.9% 


45.4% 


62.8% 


72.5% 


78.8% 


sphere 40 pc 


18.9% 


33.9% 


45.4% 


54.5% 


31.4% 


41.5% 


45.8% 


48.0% 


45.1% 


62.2% 


71.7% 


77.7% 


shell 10 pc 


9.10% 


17.5% 


25.2% 


32.2% 


41.3% 


56.9% 


65.6% 


71.3% 


45.5% 


63.0% 


72.7% 


79.0% 


shell 20 pc 


18.1% 


32.4% 


43.7% 


52.7% 


37.8% 


51.6% 


58.9% 


63.5% 


45.4% 


62.8% 


72.4% 


78.7% 


shell 40 pc 


28.2% 


45.8% 


57.0% 


64.7% 


32.5% 


43.0% 


47.5% 


49.7% 


44.9% 


61.9% 


71.2% 


77.1% 



TABLE 2 

Reference abundances 



Element 


WNM 


WNM2 


ISM 


Orion 


Bstar 


Solar 


H 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


He 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


C 


8.15 


8.15 


8.40 


8.48 


8.27 


8.39 


N 


7.88 


7.88 


7.90 


7.85 


7.80 


7.93 


O 


8.50 


8.50 


8.50 


8.60 


8.56 


8.69 


Ne 


8.07 


8.07 


8.09 


7.78 


8.20 


8.00 


S 


7.07 


7.27 


7.51 


7.00 


6.97 


7.27 



TABLE 3 

Spectral types, temperatures T» , radii R* , hydrogen ionizing luminosities Q(H°), helium ionizing luminosities Q(He°), and 

GRAVITIES g t OF IONIZING STARS 



Sp. Type 


T, (K) 


R* (cm) 


G(H°)(s- 1 ) 


G(He°) (s- 1 ) 




03V 


48,410 


1.24xl0 12 


9.73 xlO 49 


2.35 xlO 49 


3.98 


04V 


45,180 


1.13X10 12 


5.46xl0 49 


1.14xl0 4? 


3.94 


05V 


42,160 


1.03X10 12 


2.93 xlO 49 


4.95 xlO 48 


3.90 


06V 


39,350 


9.39x10" 


1.62xl0 49 


2.04X10 48 


3.86 


07V 


37,150 


8.37x10" 


8.63 xlO 48 


7.32 XlO 47 


3.85 


08V 


35,480 


7.29x10" 


4.30 xlO 48 


1.92 xlO 47 


3.87 


09V 


33,490 


5.92x10" 


1.48X10 48 


2.16X10 46 


3.95 


B0V 


31,620 


5.04x10" 


4.69 xlO 47 


2.52X10 45 


4.00 


B1V 


26,600 


4.10x10" 


1.77xl0 46 


6.70X10 42 


4.00 
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FIG. 1 . — Dust-scattered halos for various sizes and shapes of the central source. In the sphere models shown in (a) and (c), Ha photons are uniformly emitted 
from spheres with radii of 10, 20, and 40 pc. In the shell models shown in (b) and (d), Hct photons are uniformly emitted from shells with a radial thickness of 
5 pc and outer radii of 10, 20, and 40 pc. Surface brightnesses are normalized to the same total luminosity in (a) and (b), and normalized to the central surface 
brightness in (c) and (d). Black, red, and green curves correspond to sources with outer radii of 10, 20, and 40 pc, respectively. Solid and dashed lines denote 
total (direct + scattered) and scattered surface brightnesses, respectively. 
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FIG. 2. — Brightness profiles of dust-scattered halos for a spherical source with a radius of 10 pc. Hydrogen density varies from 20 cm 3 to 2 cm 3 from the 
upper to lower curves. The ordinate is shown in natural logarithmic scale to best represent the slope of exponentially decreasing function. Mean free path / ra f p 
and characteristic length r c , which was obtained by fitting the brightness profile to a function r~ 2 exp(— rj r c ), for each hydrogen density are shown in parentheses. 




Radius (pc) Radius (pc) Radius (pc) Radius (pc) 

FIG. 3. — Surface brightness profiles of OB associations in R-band image. Model profiles for OB association with radii of (a) 10 pc and (b) 30 pc, which are 
expected from a telescope with a spatial resolution of 1.2". Black curves simulate brightness profiles of OB associations in the presence of diffuse background. 
The expected profiles in the case of no background are shown in red color. Hydrogen density was 10 cm -3 for the models. Four sample profiles of OB associations 
in NGC 7793, taken from SINGS data, are shown in (c) and (d). In the figures (c) and (d), the coordinates (Aa, A<5) of two OB associations relative to the galactic 
center (a, 8) = (23 h 57 m 49.83 s , -32 d 35 m 27.7 s ) are shown in units of arcmin in parentheses. 
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Radius (pc) Radius (pc) Radius (pc) 

FIG. 4. — Top: brightness profiles of Ha, He I A5867, [N II] A6583, and [S II] A6716 lines (in units of erg cm 2 s~' sr -1 ) for various central ionization sources. 
Bottom: brightness profiles of line ratios He I/Ha, [N II]/Ha, and [S II]/Ha. Elemental abundances for WNM and hydrogen density of nu = 10 cm~ 3 were 
assumed for the photoionization models. The curves from the outermost to innermost correspond to 03V to B1V stars progressively. Solid and dashed lines were 
alternatively used for clarification. 
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FIG. 5. — Brightness profiles of dust-scattered optical emission lines originating from the H II region surrounding 04V (solid lines) and 09V (dotted lines) 



stars. Top panels from left to right show intensity profiles of Ha, He I A5867, [N II] A6583, and [S II] A6716 lines (in units of erg cm 



sr ). The intensities 



for an 09V star are scaled up by a factor of 2. Bottom panels from left to right show profiles of line ratios He I/Ha, [N II]/Ha, and [S II]/Ha. Red and blue lines 
represent direct and scattered emission lines, respectively. Thick black lines are total (direct + scattered) intensities or intensity ratios. Elemental abundances for 
WNM and hydrogen density of 10 cm -3 were assumed for the calculation of the H II region. The extent of the dust cloud in the z-direction was \z\ < 100 pc. 
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FIG. 6. — Luminosity ratios of He I/Ha, [N II]/Ha, [S II]/Ha, and [S II]/[N II] for various central ionization sources (from 03V to B1V stars), abundances 
(WNM, WNM2, B star, ISM, and Orion nebula) and hydrogen densities (0.1, 1,0, 10, and 100 cm -3 ). The photoionization models were calculated with CLOUDY. 
Solid and dotted lines were obtained by using the Kurucz and WMBASIC atomospheric models, respectively. 
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FIG. 7. — Luminosity ratios of He I/Ha, [N II]/Ha, [S II]/Ha, and [S II]/[N II] for various central ionization sources (from 03V to B1V stars), abundances 
(WNM, WNM2, B star, ISM, and Orion nebula). Hydrogen density is 10 cm~ 3 . Top, middle, and bottom panels represent the line ratios obtained from the models 
"U," "CI," and "C2," respectively. The models were calculated with MOCASSIN. 
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FIG. 8. — Line ratios ([N II]/Ha, and [S II]/Ha) in the dust-scattered halos of H II regions surrounding 04V and 09V stars. Average hydrogen density of 10 
cm -3 and the WNM abundances were assumed. The left four figures represent line ratios of [N II]/Hq, and the right four figures [S II]/Ha. Numbers indicate the 
intensity ratios at the dust-scattered halo outside of H II regions. 
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FIG. 9. — Comparison between the line ratio maps (He I/Ho, [N II]/Ha, and [S II]/Ha from left to right) obtained from the uniform (top) and clumpy (bottom) 
models. 
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FlG. 10. — Increasing factor of the line ratios [N II]/Ho and [S II]/Hq, due to the underlying H« absorption, as a function of distance from an H II region. Here, 
r o(~ Wp) i s the distance from which Equation 1 is applicable. EW of Hct absorption line relative to the diffuse continuum and EW of Hct emission line relative 
to continuum of OB association are denoted in parenthesis , WS^,). 



